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^ ■ Abstract 



A new resonance version of NLS equation is found and embedded to the reaction- 
diffusion system, equivalent to the anti-de Sitter valued Heisenberg model, realizing a 
particular gauge fixing condition of the Jackiw-Teitelboim gravity. The space-time points 
where dispersion change the sign correspond to the event horizon, and the soliton solutions 
to the AdS black holes. The soliton with velocity bounded above describes evolution on 



(N 

. the hyperboloid with nontrivial winding number and create under collisions the resonance 

^ ! states with a specihc life time. 

m 
o 

00 ' !• Resonance NLS equation. 



The connection between black hole physics and the theory of supersonic acoustic flow 
was established by Unruh 1 and has been developed to investigate the Hawking radiation 



and other phenomena for understanding quantum gravity . Recently, the similar idea for 
simulation of quantum effects related to event horizons and ergoregions but in superflu- 
ids, which in contrast to the usual liquids allow nondissipative motion of the flow, was 
^ proposed 3 . In this case a " superluminally" moving inhomogeneity of the order param- 

eter like solitons, provides black holes like quasi-equilibrium states, exhibiting an event 
horizon. Although Jacobson and Volovik considered a simplified profile of soliton and 
mentioned unimportance of exact structure of the solution, an exactly soluble model of 
solitons related to black holes allows one to describe the scattering process of black holes 
and corresponding quantum phenomena. Very recently, some attempt was done to describe 
solitons of integrable models like the Liouville 4 , the Sine-Gordon 5 ' 15 and the Reaction- 
Diffusion system 6,7 , as black holes of Jackiw-Teitelboim (JT) gravity. In the last paper 
the scattering of two identical soliton-like structures called dissipatons, relating to the 
black holes and creating the metastable state was considered. Continuing in this direction, 
in the present paper we present a new integrable version of the Nonlinear Schrodinger 
equation (NLS) in 1+1 dimensions, admitting solitons with a rich resonance scattering 
phenomenology and interpreted as black holes of JT gravity. 
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The NLS, well known also as the Gross-Pitaevskiy equation, appears in the phe- 
nomeno logical description of superfluidity of an almost ideal Bose gas 8 . In this case, the 
squared modulus of the wave function tfjtfj is interpreted as the particle number density 
in the condensate state, while the gradient of the phase is proportional to the superfluid 
velocity v = Varg-i/;. The similar hydrodynamical interpretation for the decomposition of 
quantum mechanical wave function was considered by Madelung long time ago 9 . It can be 
extended to the NLS, representing a generic envelope equation and appearing in the wide 
range of phenomena, with many applications. Besides superfluidity, the most popular on 
is probably in the nonlinear optics, describing light propagation in optical fibers 9 . The de- 
composition shows that quadratic dispersion term in the NLS results from contribution of 
two terms: the superflow (classical drift) dispersion and the so called "quantum potential" 



The potential was introduced by de Broglie and explored by D. Bohm 11 to make a 
hidden-variable theory and is responsible for producing the quantum behavior, so that 
all quantum features are related to its special properties. In more recent context it is 
considered in the stochastic mechanics as producing non-classical diffusion 12 . Relation of 
non-classical motion with the internal spin motion or the zitterbewegung was established 
in a series of papers 13 . In the optical context the superflow term corresponds to the 
geometrical optics, while the quantum potential to the diffraction, which leads to the 
spontaneous pattern formation and allows the quantum effects even at room temperature 
and macroscopic scale. We emphasize that existence of the envelope solitons for nonlinear 
equations is possible only with this potential, since its dispersive contribution is exactly 
compensated by the nonlinearity. Changing the sign of quantum potential contribution to 
the NLS, modify the behavior of solitons drastically. The potential Uq is invariant under 
the complex constant rescaling transformations, ip(x,t) — > Xt^(x,t), A G C, and thus does 
not depend on the strength of the wave, associated with soliton , but depends only of the 
form of the wave. Therefore, its effect could be large even for the well separated and far 
enough solitons. 

Below we consider the NLS soliton propagation in the quantum potential (1.1): 



which we call the resonance Nonlinear Schrddinger equation (RNLS). It could corresponds 
to the response of a hypothetical resonance medium to an action of a quasimonochromatic 
wave with complex amplitude ip(x, t), which is slowly varying function of the coordinate and 
the time. Eq.(1.2) can be considered as the third integrable version of NLS, intermediating 
between the defocusing and focusing cases (repulsive and attractive non-ideal Bose gas 
correspondingly). The additional term in the right hand side of NLS (1.2), depending 
on the form of the amplitude profile of the wave, can be interpreted also as a specific 
electrostriction pressure or the diffraction "force" 14 . Worth to note that the model (1.2) 
is integrable only with the chosen coefficient 2 on the r.h.s, but the resonance properties 
can be valid as well in a more general case. 




(1.1) 




(1.2) 
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2. Madelung fluid and the reaction-diffusion system. 



By decomposing the wave function ip e C , ip = ex.p(R — iS), ip = exp(i? + iS), in 
terms of two real functions R, S G 1Z, the model (1.2) can be represented as the Madelung 
fluid 9 

-d t R + d 2 x S + 2d x R d x S = 0, (2.1a) 
-d t S + d 2 x R + (d x R) 2 + (d x S) 2 - je 2R = 0. (2.16) 
Introducing two new real functions e + = exp(i? + S), — e~ = exp(R — S 1 ), or 

_ e+e - = e » = M ., s= i ln ^ = i. ln | (2 . 2) 

we have the system 

-d t e + + d 2 e + + ^e + e~e + = 0, (2.3a) 

+d t e~ + d 2 e~ + je+e'e' = 0, (2.36) 

representing a particular form of a 2-component reaction - diffusion (RD) system. We 
note that unusual negative value for diffusion coefficient in the second equation is crucial 
for the existence of Hamiltonian structure and the integrability of the model. In this 
case the system (2.3) is time reversible t — > — t and invariant under the global SO(l, 1) 
transformations — > e ±a e ± . It admits solution with exponentially growing and decaying 
components, 

e ± = ±( A)^ e ±[(^ 2 +fc 2 )i-|H cosh" 1 ^ -vt- x )}, (2.4) 
but with perfect solitonic shape for the 0(1, 1) scalar product 

-e + e~ = |V>| 2 = —^k 2 cosh~ 2 [k(x - vt - x )}. (2.5) 

By analogy with the dissipative structures in the pattern formation we called these dissi- 
pative soliton solutions as dissipatons 6 . 

Using (2.2) we can find one-soliton solution of Eq.(1.2) corresponding to the one- 
dissipaton solution (2.4) 

^=(^-)2 k - — — . (2.6) 

As we see the role of the "quantum potential" is to change the dispersion for the NLS. 
Indeed, when Uq = 0, Eq.(1.2) is the so-called defocusing NLS, admitting the "dark" 
soliton solution with non-vanishing boundary values. While with Uq ^ 0, we have the 
"bright" soliton (2.6), with vanishing boundaries, which usually appears in the focusing 
NLS. 
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3. The gravitational interpretation. 



Dissipatons are related to the black hole solutions of the Jackiw-Teitelboim gravity 4 ' 7 
and provide interesting tools to study nonperturbative sector of the General Relativity. 
Defining the 2-dimensional metric tensor in terms of Einstein-Cartan zweibeins 



9i*» = efelVab = ^(e+e u + e+e^ ), (3.1) 
where e J = ej ± e 1 = (ej, ef ), i] ab = diag(-l, 1), and 

4 = ±^e±, ef = e ± , (3.2) 



such that 



de+ de + _ 1 de+ _ +de 

900 = ~^^ 911 = 6 6 ' 901 =2^ ~ e 0^' (3 - 3) 



we find that for e^ 1 satisfying (2.3) the metric describes two dimensional pseudo - Rieman- 
nian space-time with constant curvature A: 

R = g^Rw = A- (3.4) 

This, low dimensional (lineal) gravity model is known as the Jackiw-Teitelboim model 4 . It 
can be used also to describe the S-wave sector of the extremal D = 4 supersymmetric black 
hole solutions of models with specific dilaton coupling 15 . When the curvature vanishes, 
A = 0, the nonlinear term in Eqs.(1.2) and (2.3) disappears. Then the system (2.3) 
becomes the decoupled linear heat equation and the time-reversal one, while (1.2) reduces 
to the unusual, sign indefinite dispersive modification of the linear Schrodinger equation: 

i^ + fl^-2^V = 0. (3.5) 

Worth to note that the classical theories corresponding to the usual quantum mechani- 
cal Schrodinger equation and to our modified model (3.5) would be equivalent, since the 
"quantum potential" is proportional to h, and in the % — > limit both models lead to the 
same Hamilton- Jacobi equations. It turns out that this model is relevant to the black hole 
solutions of the CGHS string-inspired gravitational theory 16 . 

The gauge fixing condition (3.2) is equivalent to the classical -SO (2, 1) / -SO(l, 1) Heisen- 
berg model on the anti-de Sitter space (A < 0), 

d Q s = s A d 2 x s, (3.6) 

where are local coordinates in the tangent plane <9 M s = (^)^(e+n_ — e~n + ), so that 
the metric tensor is 

2 

9»v = (—r)(d»sd v s), (3.7) 
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and s 2 = -(S 1 ) 2 + (S 2 ) 2 - (S 3 ) 3 = -1. Actually in Sec.7 we show that (3.6) is gauge 
equivalent to (2.3). 

Formulas (3.3) allow us to establish a correspondence between the geometrical and 
physical characteristics of the model. Indeed, in terms of ip variables the metric tensor is 
given by 

so that goo component is the dispersive part of energy density, while gn and goi the mass 
and momentum densities correspondingly. The mass, momentum and energy conserved 
quantities 

/oo p oo 

e + e~dx= I \ip\ 2 dx, (3.9a) 
-oo J — oo 

/OO />CO 
(e + d x e~ — d x e + e~) dx = i / (d x i(ji(j — ^d x i\i) dx, (3.96) 
-oo J —oo 



pOO A pOO _ A 

E = 2 / [d x e+d x e---(e+e-) 2 }dx = 2 / [d x ^d x ^-2d x \^\dM - -\^\ 4 ] dx, (3.9c) 

J-oo O J-oo O 



for one-dissipaton solution (2.4) or the one-soliton solution (2.6) are 

16 Mv 2 A 2 

M=—\k\, P = Mv, E = — + —M 3 . (3.10) 

These expressions show that it can be interpreted as non-relativistic quasi-particle of non- 
negative mass M and momentum P, with the positive rest energy Eq = E(v = 0) = ^M 3 . 

4. Resonance dispersion and black holes. 

As we will see below an interaction of these particles leads to the creation and annihi- 
lation processes, forming resonance states. The existence of these states follows from the 
form of dispersion part of energy density (3.9c), written in terms of variables (2.2) 

e = 2{d x i>d x ^ - 2c\#|a i #|) = 2[{d x S) 2 - (d x R) 2 ]e 2R . (4.1) 

For the NLS case, when Uq = 0, we have the positive definite dispersion energy 

e = 2d x 4>d x ^ = 2[{d x S) 2 + (d x R) 2 ]e 2R , (4.2) 

while in contrast for Uq ^ 0, from (4.1) it becomes indefinite. Function eo changes the 
sign at space-time points where 

(d x S) 2 - (d x R) 2 = (d x S - d x R)(d x S + d x R) = 0, (4.3) 

or d x S = rkd x R. Now, if we compare (4.1) with (3.8), so that the dispersion energy density 
has geometrical meaning of the metric tensor component eo = — 2(7oo, then conditions (4.3) 



are equivalent to existence of the event horizon space-time points (##,£#), where goo, 
change the sign. In this way we relate the resonance dispersion of NLS with the " quantum 
potential" (1.2), to the event horizon in two dimensional space-time. The existence of 
event horizon indicates the nontrivial causal structure of the space-time and the black hole 
phenomena. In fact, if we calculate the metric (3.8) for one-soliton solution (2.6) 

ds 2 = -^-[(A; 2 tanh 2 k(x - vt) - ^v 2 )(dt) 2 - (dx) 2 - vdxdt]\tjj\ 2 , (4.4) 

it shows a horizon singularity at 

v 

tanh/c(x - vt) = ± — , only if, \v \< 2 \k \ = \v max \. (4.5) 
2k 

Consequently, a black hole dissipaton cannot move faster than the maximal value of the 
velocity liWxl = 2|fc|. In this case the metric can be transformed to the Schwarzschild 
type form and shows the causal structure in terms of Kruskal-Szekeres coordinates 7 . 



5. Hydrodynamical interpretation. 



The Madelung fluid representation gives simple hydrodynamical explanation of of the 
resonance states and the event horizon existence. In (2.2) representation we introduce the 
fluid density p = e 2R = \i(j\ 2 and the local velocity V = —2d x S, describing the center 
of mass motion. Then, characterizing the internal motion referred to the center of mass 
frame, the stochastic diffusion or the zitterbewegung 13 non-classical contribution from the 
"quantum potential", can be described by velocity Vq = For the normal Madelung 
fluid (the NLS, when Uq = 0) the dispersion energy density (4.2) is just the sum of kinetic 
energies of these two motions 

eo = (-^- + -^)- (5.1) 
In contrast, for the non- vanishing Uq ^ 0, the density (4.1) is the difference 

pV 2 pVZ 

eo = ( ^r- 2 } - (5 - 2) 

In this hydrodynamical representation, the metric tensor (3.8) has the form 

9oo = \p( v q - y2 )i 9ii = ~P, 9oi = ^pV, (5.3) 

which is similar to the ADM split of a (l+l)-dimensional Lorentzian spacetime correspond- 
ing to the so-called acoustic metric, derived by Unruh 1 for the sound waves in a fluid 2 . 
Then, the existence of resonance states or the event horizon has the meaning of the equality 
of the center of mass and internal motion velocities V = ±Vq. For the one-soliton solution 
(2.6) velocity of the classical center of mass motion is constant and coincides with the soli- 
ton propagation velocity V = v, while the " quantum" velocity is Vq = —2k tanh k(x — vt) . 
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The last one is bounden above by the constant value \Vq\ < 2\k\, why the velocities com- 
pensation is possible only in this region. Then, the compensation condition V = ±Vq is 
equivalent to the event horizon Eq. (4.5). 

To derive the black hole we first rewrite the metric (5.3) in the moving frame (£, t) = 
(x — vt, t) with constant velocity v. Then it has convenient form in terms of the new 
"shifted" local velocity W(£,t) = 2v - V(£,t): 

9oo = \p{Vq - W 2 ), gn = -p, g i = ~pW. (5.4) 



For a one-soliton solution (2.6) or dissipaton (2.4), moving with a constant velocity v, 
we have W = v and the metric (5.3), where V = v = const., but with stationary space-time 
geometry due to the time independence of p = p(£), Vq = Vq(£). 

In the general case the metric contains the off-diagonal terms. The time synchroniza- 
tion of this space-time is possible if function 2W/(Vq — W 2 ) is integrable. Then we can 
define new time coordinate dr = dt — 2W/ (Vq — W 2 )d^ and get the metric 



ds 2 = P [\{V^-W 2 ){drf-^^{de] 

Q 



(5.5) 



Following the same arguments as for the black holes, the Hawking temperature can be 
derived from this metric. In particular case of one soliton solution everything can be done 
explicitly. Synchronization of the stationary metric is given by the above time transfor- 
mation integrated as 



r = m,t) = t + 



2fc 3 (l- 7 2 ) 



1 



2|7l 



In 



I7I + tanh k£ 



I7I — tanh k£ 



where I7I = \v/2k\ < 1. The Hawking temperature in this case is Th = -^k 2 (l — r y 2 ). 

Comparing (5.2) with (3.9c) and (3.10) we can see that the kinetic part of the internal 
motion gives the negative sign contribution to the rest energy Eq, but due to the positive 
potential part, the resulting energy is positive. This is the reason why we can have the 
resonance behavior for our model, but not for the NLS. Indeed, decay of a soliton in the 
rest means creation of a pair with the positive energy, which is allowed only if the rest 
energy is positive. From conservation of the mass (3.9a) follows that the mass defect for 
soliton decay is always zero, AM = M — (Mi + M 2 ) = 0. Then, the rest energy satisfies 
the condition 



En = 



A 2 A 2 

M 3 = (Mi + M 2 Y 

384 384 v ; 



> 



A 2 ^ 
384 



(Mf + Ml) =E 0{1) + E m , (5.6) 



such that AE Q = Eq — (-E'o(i) +-^0(2)) > 0, which allows creating two solitons. In contrast, 
in the NLS case, the sign of the rest energy is negative and instead of inequality (5.6) one 
has Eq < -E'o(i) + Eq(2) , which forbids decay of the bright soliton. 
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6. Resonance interaction of solitons. 



Rewriting RD system (2.3) in the bilinear form 

(±D t -Dl)(G ± oF) = 0, D 2 x (FoF) = -2G+G-, (6.1) 

where three new real functions are defined by = (— 8/A)^G ± /F, with the product 
e + e~ = — zr/{(logF) xx , we apply the Hirota bilinear approach to construct solutions of 
(2.3). 

The one-dissipaton is given by following solution of the system (6.1) 

G ± =±e"i ± , F = l + e ^+%+^.i, e^ 1 - 1 = (k+ + fcf)" 2 , (6.2) 

where rjf = /c^x ± (kf) 2 t + , and in terms of redefined parameters, fe = (&/" + fcj~)/2, 
w = — (fej" — fc]~) it has the form (2.4). 

Now we compare the one-dissipaton boundary conditions with the horizon condition 
(4.5). In the space of parameters (v, k) there exist the critical value v cr u = 2k. For solution 
(6.2) when v < v cr u, one has e^ — > at infinities. So the vanishing b.c. for dissipaton are 
equivalent to the black hole (BH) existence. The corresponding heavy particle we call the 
BH dissipaton. At the critical value the solution is a kink steady state in the moving frame 
e ± = ±ke ±k ^°(l =F tanh k£), with constant asymptotics e ± — > ± 2ke ±k ^° for x — > =Foo 
and — > ± for a; — > ±oo. In this case we have the extremal black hole or the EBH 
dissipaton. In the over-critical case v > v cr i t , e ± — > ±oo for x — > =Foo and e^ — > ±0 for 
x — > ±oo, no black hole exists and we have very fast and light particles called the LD. 

For the two-dissipaton solution we have 

G ± = ±[eit + + (-^^fe^+^+vf + { ^_f e vt^^ ] . (6 . 3a) 

k 2 i k 1± k 12 k 2 2 

evt+Vi e vt+V2 e v}+vT e v}+V2 



F = l + —— + —— + 



(ktn 2 (kt 2 ~r (tin 2 (* 



+- 

22 



"a2 ^21 "'ll ^22 

where fc$ = fef + fej, k^ = kf - fcj, 7/f = fef x ± (Jfef ) 2 t + t? ±(0) . 

First we consider the degenerate case of (6.3), when kf = k^ = pi, k 2 = k 2 = P2- 
Then the solution can be simplified and after some transformation be represented in the 
following form 

_i_ ,,8.i pi cosh^ 2 e ±p i* +P2 cosh^!e ±p 2* 

6 — i( ) 2 V+T>— (6.4) 

—A p 2 _ cosh 9 + + p\ cosh 9- + 4pip 2 cosh(p + p_t) ' 

where, p± = pi ±^2, 9± = 6 1 ±6 2: 0i = pi(x — x 0i ), (i = 1,2). It describes collision of two 
identical dissipatons with amplitudes p+/2, moving with the equal velocities \v\ = \p-\ 



and creating the resonance state with the life time, AT 2p 2 d/p + p-, depending of the 
relative distance d, where xoi = 0, xq 2 = d. This particular solution was considered in 
Ref. 7. 

In a more general case, when kf > 0, (i = 1,2), and fe^T" > 0, fc^ - > 0, k± 2 ~ > 
0, k^i < 0, (6.3) describes collision of two different dissipatons with amplitudes fc^~~ /2 
and fc^ - '/2 and velocities U12 = —k^ 2 ~ and i^i = — fc^i - correspondingly. Depending of 
the positions shift the resonance states can be created. 

As a simplest case we consider conditions for decay of a BH dissipaton at the rest 
(v = 0) on two dissipatons (2.4) with parameters (ki,vi) and (k 2 ,V2)- From the mass, 
momentum and energy conservation laws we obtain relations 

v 2 = 4k 2 , v 2 = 4k 2 . 

Two possibilities exist: 

a) \vi I = \v 2 1- In this case |fci| = j/^l and both dissipatons have the equal mass 
Mi = M2 = M/2 and velocities, satisfying the extremal conditions vf = 4kf, (i = 1,2) 
and corresponding to two EBH illustrated on Fig.l (for t > 30). 

b) \vi \ > \v 2 \ (without lose of generality). In this case vf > Ak\ and v\ < 4/c|, and we 
have decay on the one BH dissipaton and the one LD as illustrated in Fig.l (for t < 0). 

On Fig. 2 the interaction of two BH dissipatons by exchange of LD is shown. A more 
complicated interaction (Feynman diagram) with 4 vertices is shown on Fig. 3. Detailed 
description of various interactions simulated by MATHEMATICA would be published 
elsewhere. 



7. Black holes as topological solitons. 



The gauge equivalence between (2.3) and (3.6) allows us to construct exact solutions of 
(3.6) providing simple geometrical visualization of the event horizon position and to treat 
the black hole as topological soliton. Under collision they should have similar resonance 
properties as solutions of Eqs. (1.2) or (2.3). In the matrix representation for 50(2, 1) 

S = i(a+ S q3 J = ( s > r ) = arsg' 1 , (7.1) 



x s+ -s\ 

where S ± = S 1 ± S 2 , S 2 = —I, detS = 1, Eq.(3.6) has the standard matrix form 

d t S=±-.[S,d 2 x S], (7.2) 

with the Lax pair 

J? M = 7 AS, J HM = l -X 2 S + jSd x S. (7.3) 

4 o 4 

This model is gauge equivalent to the resonance NLS (1.2) and RD (2.3). Although the Lax 
pair for the first one can be derived 17 , for the second one it has simple Zakharov-Shabat 
form 

T RD_f\^ q~ \ T RD_f^ 2 -q + q~ -(d x -±\)q-\ (7 
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where = (-jp) ^e^. The gauge equivalence is implemented in the standard way by non- 
Abelian transformation J^ M = gJ^ D g~ x — d^gg -1 , (// = 0, 1) where the matrix g(x,t) 
is solution of the linear problem d^g = gJ^ D (X = 0). To construct the magnetic analog 
of dissipaton solution (2.4) we solve this system first. The result is 

, , / tanh^-7 _i e 7*-fc 2 (i-7 2 )t\ 

9(x,t)=[ . ,J t , 2 ,. coshz , (7.5) 

\ l „- 7 ^+fc (1-7 )* tanh^ + 7 ' 



cosh 2 



where z = k(x — vt), 7 = v/2k, detg = 1 — 7 2 . Thus, from (7.1) we find one magnetic 
dissipaton solution with components 

S 3 = -1 + 2 _ cosh" 2 z, (7.6a) 
1 — 7 2 

5" = cosh" 1 z (tanhz - 7) e ^-^^)\ (7.66) 

5+ = — ^ cosh" 1 z (tanhz + 7) e -7*+fc 2 (i-7 2 )t. ( 7 . 6c ) 
1 — 7 2 

It describes magnetic (curved) analog of dissipaton (2.4). Indeed, in the moving frame 
with velocity v the S s component is time independent while the S~ and S + components 
are decaying and growing in time. As well as for dissipatons, properties of solution (7.6) 
essentially depend on the velocity v. We have following cases. 

a) 7 2 < 1, or \v\ < 2\k\, then -1 < S 3 < At any fixed time t, S + and S~ 
vanish when z — > ±00, while S 3 — > —1. Due to s = (0, 0, —1) when z — > ±00, the real line 
R is compactified. Since the hyperboloid s 2 = — (5" 1 ) 2 + (>S 2 ) 2 — (S 3 ) 2 = —1 has topology 
of cylinder R x S 1 , the solution (7.6) describes the mapping 5" 1 — > 5" 1 with degree one. 
Thus, we see that (7.6) is the topological soliton, travelling with constant velocity v. When 
z = zh, so that tanhz = ±7, one of the components S + or S~ is zero. According to (3.7), 
in this case the metric component goo = and like dissipatons we have the event horizon at 
(4.5) (see Fig. 4). Since any topological soliton configuration cross each of the lines S + = 
and S~ = at least once, the intersection points correspond to the event horizon. This 
fact relates the black holes solution with topological soliton on hyperboloid. 

b) 7 2 > 1, or \v\ > 2\k\. In this case < S 3 < -1. At z — > ±00 one of 
the components S + or S~ is exponentially growing. Moreover, asymptotics at +00 and 
—00 are orthogonal. This is why the solution (7.6) has trivial winding. In fact it can be 
considered as the turning travelling wave, which never cross the asymptotic lines at finite 
distance. Thus, in this case there is no event horizon and black holes. 

From the above analysis we can see that existence of black holes and event horizon 
intimately related to the topologically nontrivial solutions of the model (3.6). 

8. Conclusions 

The black hole picture with resonance interaction described above can be applied to 
physical models of slowly varying quasimonochromatic wave in nonlinear media with the 
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sign indefinite dispersion. The 1+1 dimensional case can be realized particularly in the 
nonlinear optics. The great variety of optical solitons is due to many different properties 
of the media involved, including nonlinearity, material and geometric dispersion, passive 
or active properties, etc. The crucial role in soliton properties play the group- velocity 
dispersion of the optical fibers, which depends not only on the property of glass material, 
but also on the waveguide property of the fiber. According to the sign of dispersion (positive 
or negative(anomalous)), two types of NLS are known - defocusing and focusing cases, 
which admit the "dark" and the "bright" solitons correspondingly. Since the envelope 
wave function is a complex quantity, the quadratic dispersion in general consists of two 
parts: the phase dispersion and the modulus dispersion. The first one corresponds to 
the geometrical optics, while the second one is responsible for the diffraction. But in 
both focusing and defocusing cases, contribution to dispersion from the phase and the 
modulus has the same sign (positive and negative correspondingly). In the present paper 
we considered the nonlinear media with the sign indefinite quadratic dispersion, which is 
the result of competition between contributing with the opposite signs of the phase and 
the modulus dispersions. Not going into microscopic details we confined ourself here only 
by phenomenological description of some hypothetical medium. Then, the response of the 
medium to an action of a quasimono chromatic wave with complex amplitude if)(x, t), which 
is slowly varying function of the coordinate and the time, is described by our resonance 
version of NLS. If it can be realized experimentally, the resonance collisions of solitons 
would be an interesting tool in optical communication systems with the black hole physics 
flavor. 
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Figures 



Fig. la. 3D plot of two dissipatons resonance- type collision for /c+ = 0.1, k t = 1, k£ = 1, 
k% = and d = 30 in the (x,t) plane. 



Fig. lb. Contour plot of two dissipatons collision with BH resonance in the (x,t) plane. 



Fig. 2a. 3D plot of two BH dissipatons exchange- type collision for k± = 2, k 1 = 1, k£ = — 1-7, 
fcj = —1.9 and d = 50. 



Fig. 2b. Contour plot of two BH dissipatons exchange-type collision in the (x,t) plane. 



Fig. 3. Contour plot of two-BH dissipatons 4 vertex- type collision for k± = 2, fcj = 1, k£ = 1, 
k^ = 0.3 and d = 30 in the (x,t) plane. 



Fig. 4. Parametric plot of topological soliton projection on (S , S 2 ) plane for k = 1, 7 = 0.5. 

Positions of the BH event horizon correspond to intersection points with S + = and 
S~ = lines, while asymtotics on ±00 to the begining of coordinates (0, 0). 
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